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ABSTRACT 

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- 
selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies 
were first identified with Herschel/SFIRE imaging data from the Herschel Astrophysical Terahertz 
Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The 
targets were observed with Spitzer/IRAC at 3.6 and 4.5/xm. Due to the spatial resolution of the IRAC 
observations at the level of 2", the leasing features of a background DSFG in the near-infrared are 
blended with the fiux from the foreground leasing galaxy in the IRAC imaging data. We make use 
of higher resolution Huhble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles 
of the foreground leasing galaxy (or galaxies) as a way to model the foreground components, in order 
to successfully disentangle the foreground lens and background source flux densities in the IRAC 
images. The flux density measurements at 3.6 and 4.5//m, once combined with Hubble/WFC3 and 
Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution 
of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs 
to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters 
for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 
8 X 10^*^ to 4 X 10^^ Mq and star-formation rates of around 100 M© yr“^. This puts these lensed 
sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies 
at similar redshifts. The high values of SFR inferred for these systems are consistent with a major 
merger-driven scenario for star formation. 
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1. INTRODUCTION 

Dusty star-formi ng galaxies (DSF Gs; For a re¬ 
cent review, see iCasev et al.l l2014l i are now be¬ 
lieved to significant contributors to cosmic star- 
form a tion in the early Un i verse (e.g. JL^ Floch et’all 
2005t iPerez-Goiizalez et alJ 120051: iTakeuchi et alJ 120051 : 
Marchesini et al.l 1201411 . The extreme star-bursting ex¬ 
amples of DSFG s appear as bright s ub-millimeter galax¬ 
ies (SMGs, s ee iSmail et ^ Il997t iHughes et alJ 119981 
Barger et alJ Il998t iCoppin et^TI 120081 : lAustermann 
2009H and are best studied at far-IR and sub-millimeter 
wavelengths due to the high dust extinction at rest-frame 
optical wavelengths. Such galaxies have star-formation 
rates in excess of 100 Mq yr“^ and have emissions peak¬ 
ing in the far-infrared, with luminosities Tfir > 10^^ 
Lfn. The strong clustering of these starburst galaxies 


(IBlain et al.l 

2004lFarrah et alJl2006l:ICoorav et alJl2010l: 

iHickox et al. 

I2OI2D is such that they may evolve into dark 


matter halos that host some of the most luminous and 
massive elliptical galaxies today. Thus their rapid forma¬ 
tion may have moved these galaxies to the r ed sequence 
that already exists for galaxies at z > 1 (jFaber et al.l 
120071: iBarro et~aI1120131 : lllbert et al.l 120131) . It is gener¬ 
ally believed that a large fraction of the high-redshift 
DSFGs ar e trigged by galaxy merg ers, similar to local 
ULIRGS (|Sanders fc Mirab^ 11996 ). There is also ob¬ 
servational evidence for lack of merger features in the 
optica l and infrared images i n > 40% of the DSFGs at 
z ^ 1 (|KartalteDe et al.ll201lD . 

Despite significant progress in understanding the for¬ 
mation and evolution of DSFGs at 2 > 1, it still remains 
difficult to measure basic properties, such as the total 
amount of gas, stars, and dust in these objects, as well 
as the spatial distribution of these quantities within the 
galaxies. The primary obstacles for such detailed ob¬ 
servations come from the poor spatial resolution of the 
far-infrared/sub-millimeter observations with which the 
DSFGs are identified, and their faintness at near-infrared 
wavelengths. The former makes it challenging to identify 
their counterparts at other wavelengths, while the latter 
makes follow-up observations time-consuming. Never¬ 
theless, the rest-frame near-IR flux densities of galaxies 
are crucial to establish the total stellar mass content and 
to study the spatial distribution of stellar populations in 
such objects. While these studies may be challenging for 
typical dusty starbursts, gravitationally lensed DSFGs 
provide a mechanism to partially overcome the limita¬ 
tion associated with the faintness of the near-IR counter¬ 
parts, specifically through the flux enhancement associ¬ 
ated with leasing magnification. Lensed DSFGs also pro¬ 
vide enhanced spatial resolution through lensing magni¬ 
fication, and this enhancement has provided information 
down to 200-300pc scales within the interstellar medium 
of some of the lensed SM Gs, such as S MMJ2I35-0I02 
ofiSwinbank et all (|20I0D and SDP.81 of iNegrello et al.l 
l|20IOD that has been studied with an extended observa ¬ 
tion using long baselines with ALMA Ibve et all (120151) : 
IPartnership et al.l (|20I5D . 

Despite progress in high-resolution imaging in the mm- 
wavelengths with interferometers and large optical and 
infrared telescopes in the rest-frame wavelengths, our 
views of distant galaxies at wavelengths above 2.2 /rm 
and above in the mid-infrared region are still limited to 


imaging data that are diffraction limited. For galaxies 
at z ~ 1 to 3, such data are still crucial since they pro¬ 
vide necessary information to break certain degeneracies 
in models of the spectral energy distribution (SED). A 
key ingredient from SED analysis is the stellar mass of 
the galaxies, and with data out to 2.2 /rm only, it is gen¬ 
erally hard to accurately estimate the stellar mass. One 
issue is that due to high extinction, many of the DSFGs 
are faint in the rest-frame UV observable at wavelengths 
below 1 /rm. Despite limitations, Spitzer/IRAG data re¬ 
main key to estimating the stellar masses of DSFGs, as 
past studies have demonstrat ed (iHopwood et al.l 120111 : 
liviichalowski et alll2012L I2014D . 

In this paper we discuss the rest-frame optical to sub¬ 
millimeter spectral energy distributions of six lensed DS¬ 
FGs. Our key measurements in the optical to infrared 
are with Spitzer/IRAG at z = 1 to 3, which are crucial 
to improve the accuracies of stellar mass measurements 
using models of the spectral energy distributions. Due to 
the spatial resolution of IRAC imaging at the level of 2" 
per pixel, the lensed DSFGs are blended with their re¬ 
spective foreground lensing galaxies in our IRAG images. 
The paper is organized as follows. In the next Section 
we describe the target selection and the observations we 
have acquired. In Section |3] we outline the procedure 
that was used to de-blend the background DSFGs from 
the foreground lenses in the IRAG images. In Section |4] 
we discuss the SEDs of these galaxies and present model 
estimates of stellar mass, star-formation rate, dust mass, 
and other quantities. We conclude with a discussion of 
our results in Section [S] 

Throughout this paper we assume a Ghabrier (2003) 
initial mass function (IME) with a cutoff below 0.1 and 
above 100 Mq. We assume a AGDM cosmology with 

Hq = 70 km s“^ Mpc ^, flm = 0.3, and Ha = 0.7. 

2. TARGET SELECTION AND OBSERVATIONS 


The sample of lensed DSFGs studied in this pa¬ 
per was originally selected from two key surveys. H- 
ATLAS (lEales et al.l 120101) and HerMES (lOliver et al.l 
I2012D which were completed with the Herschel Space 
Observatory (|Pilbratt et al.l[20T^ . The lensing selec¬ 
tion at sub-millimeter wavelengths simply involves a 
flux density cut at 500p,m and an accounting for bright 
nearby galaxies or ga l axies that harbor rad io-loud AGN 
(|Negrelio et al.l 120101 IWardlow et al.l I2013D . The Her¬ 
MES lensing sample (see Riechers et ah, in prep, for 
a study on GO redshifts) and its selection is discussed in 
IWardlow et al.l (I2013D. while th e H-ATLAS sampl e is dis ¬ 
cuss ed in Negrello et al.l ^OlOD , iBussmann et al.l (|2013D , 
and iGalanog et al.l (I2014D . 

With our available Spitzer/IRAG data of several DS- 
EGs, we select those that were cl assified as “Grade A” 
lenses following the designation by IGalanog et al.l (120141 ) 
as our sources for study. These sources have visually 
obvious lensing features as observed in near-IR data, ei¬ 
ther with HST/WFGH or Keck/NIRG2 Adaptive Optics 
imaging in the 2.2/im band. Lensing morphologies 
include rings, arcs, and counter-images detected at high- 
significance. Some sources (HELS08 and HLockl2 in 
this study) are also classified as Grade A if a possible 
counter-image can be properly modeled after subtract¬ 
ing the foreground lens in the higher resolution data. As 
an extra check, DSFGs are classified as Grade A if their 
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Figure 1. 10" X 10" high-resolution data of selected DSFGs 

through various filters. Each image is oriented north upwards and 
east to the left, with tickmarks representing l". 

near-IR morphologies resemble the con figurations found 
in hig h-resolution sub-millimeter data (iBussmann et all 
Because this designation is based off detections 
in the near-infrared, there is a selection bias towards 
higher stellar masses. Our sample of six DSFGs includes 
HATL ASJ114638.0-001132 (G12v2.30) from iFu et al.l 
and HATLA SJ14 2935.3-002836 ('G15v2. 19) from 
IMessias et al.1 (|2014ll and iTimmons et ^ (1201511 . 

2.1. HST: WFC3 

Herschel-lensing candidates in the H-ATLAS and Her- 
MES fields were observed as part of the HST WFG3 Gy- 
cle 19 program (P.I. M. Negrello). Images were taken 
with the FI low filter (Ac = 1.15/im), with integration 
times typically ^ 4 minutes for each target, with a depth 
of 25.4 AB magnitude. We also include additional imag¬ 
ing observed through the F105W (Ac = l.OG^m) and 
F160W (Ac = 1.54^m) filters for HATLASJ142935.3- 
002836 that were carried out as part of a grism ob- 
servation of th a t gala xy, with grism data reported in 
I Timmons et al.l (|2015[ 1 . . g6T/WFC3 da t a red uction pro¬ 
cedures are described in IGalanog et al.l (|2014ll . 

2.2. Keck: NIRC2 LGS-AO 

Keck imaging data in the ATg-band was obtained us¬ 
ing the NIRC2 LGS-AO system during the years 2011 
to 2013. Due to LGS-AO observations, these imaging 
data have spatial resolutions comparable to or better 
than HST/WFC3 imaging data, reaching AO-corrected 
resolutions of 0.1" in the best conditions. Exposure 
times were usually 45 minutes to obtain a 5a point 
source dept h of 25.7 AB magnit ude using a 0.1" aper¬ 
ture radius (jGalanog et al.l [20l3l . Data reduction was 
carried ou t with an internal pip eline code that is de¬ 
scribed in iFu et al.l ()2012l 1201311 . Exposure times for 
high-resolution data are presented in Table [T] 

2.3. Spitzer IRAC 

The six lensed DSFGs was observed with 
Spitzer/IHKG in Gycle 8 through IRAC channels 1 
(Ac = 3.6/rm) and 2 (Ac = 4.5/xm) as part of program 


80156 (P.I. A. Cooray). Most targets were observed 
through 72 individual tiles using a 12-point box dithering 
pattern with 30 seconds of exposure each, for a total 
exposure time of 36 minutes. The exception is IHerMES 
S250 J142825.7-1-345547, which was observed through 
48 tiles for 24 minutes in total. We used the MOsaicker 
and Point source Extractor (MOPEX, vl8.5.0) program 
on basic calibrated data to carry out image reduction 
and mosaicking. Back ground estimations were carried 
out using SExtractor (jBertin fc Arnout^ 119961 1 . Final 
images were re-sampled to a pixel size of 0.6" x 0.6", 
with resolutions of 1.65" in FWHM for IRAC 3.6/im 
and 1.80" in FWHM for IRAC d.S/xm. 

3. DEBLENDING SPITZER IRAG IMAGES 

Due to 2" or worse spatial resolution provided by IRAC 
imaging, the observed DSFGs remain mostly unresolved 
in comparison to the high-resolution data (FWHM value 
of ~ 1.8" for IRAC imaging compared to ~ 0.2" for our 
HST and Keck data). This makes modeling the fore¬ 
ground and background sources difficult, due to signif¬ 
icant degeneracies between the flux distribution of the 
two components. Two previous studies on lensed DSFGs 
observed with Spitzer/IFikG de-blended the background 
source flux from that of the f oreground lens in different 
ways . iHopwood et al.1 (l20IIh used galfit (iPeng et al.l 
1201(1) directly on IRAC images to model the foreground 
lens. However, attempts to implement this procedure on 
our images were largely unsuccessful; the models were 
highly degenerate, primarily because the effective radius 
of the foreground components (~ 1-2 pixels or ~ 2-4”) 
often overlapped with the background source. The sec¬ 
ond approach—and the method that this study uses— 
is to convolve galfit models of high-resolution near-IR 
data with the IRAC beam, rescale the model into the 
same pixel scale as t he IRAC image, and subtract it 
from the IRAC data (iBothwell et al.1 (20131) . We make 
use of the higher resolution HST and/or Keck/NIRC2 
Adaptive Optics imaging to de-blend the lens flux den¬ 
sity from that of the background DSFG in our IRAC 
data. The fundamental assumption of this method is 
that these models retain the same morphologies across 
near-IR wavelengths from I.l/rm to 4.5/rm. For galaxies 
at z ~ 2 to 3 this assumption could have some issues 
since the I.l/rm images are below the 4000A break while 
3.6/im image is sampling redward of the break in the rest- 
frame optical. This is less of an issue once the ATg-band 
images are combined with IRAC imaging, since in these 
cases all of the images sample the galaxy redward of the 
4000A break. We have two targets without ground-based 
ATg-band imaging using the Keck/LGS-AO system, due 
to the lack of suitable AO guide stars within 45-60" of 
their location. In those two cases, unfortunately we are 
forced to use the HST/WFC3 I.l/rm image as that pro¬ 
vides the highest resolution imaging of the system. 

As the Keck and HST/WFC 3 data of all our tar gets 
were introduced and studied in iCalanog et al.l (I20I4D we 
used high-resolution models developed there for the de¬ 
blending process. Best-fit foreground models were ob¬ 
tained from galfit, while models for the backgroun d 
sources were found using gravlens (|Keetonl 1200II) . 
HATLASJII4638.0-001132 has been studied extensively 
in previous wo rks, and we used models provided by 
iFu et al.1 (I2012D . An empirical PSF from each IRAC im- 
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Table 1 

Summary of High Resolution Data 


lAU Name 

Short Name 

•^^source 

Ref. 

Exp. Time 

1 2 

Filter — tint ^ -Nframps 

IHerMES S250 J142825.7-1-345547 

HBootes02 

2.804 

R14 

-^Fiiow = 62 X 4, iif = 120 X 28, Ks = 80 X 27 

IHerMES S250 J171544.9-I-601239 

HFLS08 

2.264 

R14 

■Ipiiow = 62 X 4 

HATLASJ085358.9-I-015537 

G09vl.40 

2.091 

L14 

As = 80 X 45 

HATLASJ114638.0-001132 

G12v2.30 

3.259 

H14 

As = 80 X 42 

HATLASJ142935.3-002836 

G15v2.19 

1.026 

M14 

"^Fiosw = 88 X 4, Jfi60W = 62 X 4, = 120 x 10, Ks = 80 x 15 

IHerMES S250 J110016.3-I-571736 

HLockl2 

1.651 

R14 

•Ipiiow = 62 X 4 


Not e. — Redshift key: R14 = Riechers et al. (in prep); L14 = R. E. Lupu et al. (in prep); H14 = A. I. Harris et al. (in prep); and 
M14 = IMessias et al.) l|2014l 'l. Filters are Jfiosw = HS'T F105W, Jfiiow = HS'T FllOW, J^Fieow = HST FIQOW, H = Keck H-band, and 
Ks = Keck A's-band. 


tint is the exposure time per frame in seconds 

2 

A-frames is the number of independent frames 


age was obtained by stacking ~ 10 bright and unsat¬ 
urated sources in the same image, and then taking the 
average. Since the IRAC PSF is triangular, we make sure 
the PSF and the IRAC images are oriented the same be¬ 
fore convolution with the best-fit models. We normalized 
the PSF by the ratio of the peak high- to low-resolution 
fluxes. The beam-convolved near-IR models were then 
rescaled to the same pixel scale as the IRAC data and 
subtracted from the IRAC images. 

In order to obtain optimal fits and subtractions, we 
scaled the foreground and background models simulta¬ 
neously and corrected for positional errors (i.e., we fit 
for separate scale factors for both models, as well as po¬ 
sitional shifts in the horizontal and vertical directions) 
using the iterative IDL routine amoeba_sa, which uti¬ 
lizes the downhill simplex method in multi-dimensions 
to minimize multi-dimensional functions. During the fit¬ 
ting, we allow for a maximum positional shift of 1 pixel 
to prevent accidental fitting of tertiary sources, an issue 
that is especially prevalent in crowded fields, as found 
for HFLS08 and G12v2.30. Uncertainties in the scale fac¬ 
tors and position shifts were calculated using the Markov 
Chain Monte Carlo technique (MCMC). The convolved 
models and residuals are shown in the Appendix. 

We determine the goodness of the fits based on the 
minimized, reduced xi statistic, as given by: 

2 _ 1 {fdata{x,y) — fmodeljx, y))^ , , 

cr(x,y)2 

where /model is the sum of the convolved foreground and 
background models, as well as any secondary components 
(used on a case-by-case basis, see Section lHlTl) . and Fdof 
is the number of degrees of freedom (number of pixels — 
number of free parameters). 

3.1. Notes on Individual Sources 

We now briefly discuss each lensed DSFG as well as 
their individual de-blending solutions. These results are 
shown in the Appendix. 

IHerMES S250 J142825.7-|-345547 (HBo6tes02): 

The lensing morphology in high-resolution imaging 
shows an Einstein ring lensed by an edge-on foreground 
galaxy. A well-isolated source in IRAC imaging, de¬ 
blending the foreground and background lenses was 


straightforward. Light profile fits to the lensing galaxy 
are the most robust and reliable of the sample in this 
study, with xH values on the order of unity. The 
fitting at 4.5/rm was significantly worse than in 3.6/xm, 
most likely due to poor imaging and consequentially a 
suboptimal convolution with the PSF. 

IHerMES S250 J171544.9+601239 (HFLS08): 

High-resolution near-IR imaging shows an arc about 
3" east of the foreground galaxy. Moreover, this target 
is located in a crowded region, with two other sources 
a few arcseconds away to the east, whose light could 
potentially pollute the contributions from the foreground 
and background galaxies. Attempts to fit for these sec¬ 
ondary sources only introduced additional uncertainties. 
Final convolved models are highly degenerate and the 
photometry and subsequent analyses of this source 
should be treated with caution. 

HATLASJ085358.9-1-015537 (G09vl.40): The 

lensing feature is an Einstein ring centered almost 
exactly on the foreground lens. The de-blending process 
for this source is unique, because a separate fit was 
made for a source located about 2" to the southwest. 
Fitting for this object was deemed necessary because 
imperfections in the PSF introduced degeneracies and 
an increased systematic uncertainty to the light profile 
modeling. 

HATLASJ114638.0-001132 (G12v2.30): High- 

resolution data has shown the lensing morphology to be 
highly complex and lensed by four foregroun d sources 
in the near-IR—for a detailed analysis, see iFu et al.1 
(|2012ll . Similar to HFLS08, this field is fairly crowded, 
with multiple sources within 4" to the southeast and 
southwest. Attempts to fit for the nearby sources cre¬ 
ated very degenerate results, so only the foreground and 
background components were used in the de-blending 
process. It should be noted that light from neighboring 
sources may skew photometry of the background DSFG. 

HATLASJ142935.3-002836 (G15v2.19): The 

lensing morphology is an incomplete Einstein ring 
around a star-forming, edge -on spiral galaxy an d has 
been studied e x tensiv ely in iMessias et al.l (l2014fl and 
iTimmons et al.l (j2015fl . The source is sufficiently iso¬ 
lated so that no light pollution from nearby sources can 
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Table 2 


Magnification Factors of SMGs 


Mnir 

Ref. 

Mfir 

Ref. 

HBootes02 

c 0 + 1-4 
^•'^-0.4 

C14 

o 

CO 

1 + 

B13 

HFLS08 

7 7+1-6 

' •' -0.7 

C14 

8.2^ 

— 

G09vl.40 

11.41+ 

C14 

15.31+ 

B13 

G12v2.30 

00 CO 

6 6 

+ 1 

t- 

CD 

F12 

7-6li.t 

F12 

G15v2.19 

9'6la3 

C14 

1 + 

o o 

+ + 

M14 

HLockl2 

6 d 

+ 1 

o 

+ 

C14 

5.5^ 

— 


Note. — Magnification factors used to convert from observed 
to intrinsic flux densities. We use /^nir factors to convert 1.1/im 
to 4.5/^m flux densities and /ipiR factor s to convert mi l limete r flux 
densities. The reference key is: C14 = ICalanog et al.l 1120141 ). F12 
= iFiT e^^ I 112 013), B13 = IBussmann et alT1l20T^) 7'^ and M14 = 
IMessias et al'l (I20i4l ~). 

A study by ICalanog et al.l II2014I1 suggests that sub-millimeter 
fluxes are magnified by at least a factor of 1.5 greater than the 
near-lR factors, so we make this assumption here 

affect the fitting. Despite this, the final results were 
not very robust. This raises the possibility that the 
leasing morphologies observed at 3.6/tm and 4.5/im vary 
significantly from those observed in the ~ 1-2 /tm range. 

IHerMES S250 J110016.3+571736 (HLockl2): 

High-resolution near-IR imaging shows that the leasing 
morphology is an arc approximately northwest of 
the foreground source. Subtraction of the foreground 
galaxy in high-resolution imaging has also revealed 
a counter-image less than 1" to the southeast of the 
foreground lens (|Calanog et al.l l20Tll . The system lies 
in a moderately crowded field, so separate fits were 
made for the sources 4" to the southwest, 4" to the 
northwest, and 1" to the northeast. Despite the number 
of sources to fit for, the final model fit and de-blending 
was found to be robust. 

Throughout the rest of this paper, we shall refer to 
each galaxy by their respective short names for brevity. 

3.2. Photometry 

We used our own IDL codes to perform aperture pho¬ 
tometry on the IRAC background models. Data were 
converted from [MJy sr“^] to [/fJy pix“^] prior to pho¬ 
tometric measurements. We calculated the flux density 
within a circular aperture of a specified radius that would 
enclose the DSFG emissions—2" for all sources, except 
for HFLS08, where we used 4", and for G12v2.30 and 
HLockl2, where we used a 3" aperture. To find the flux 
errors, we calculated the flux densities within apertures 
of the same radii placed at empty areas of the resid¬ 
ual image. Finally, we divide the observed flux densities 
for each source by their respective magnification factors 
to find their intrinsic flux densities. We assume con¬ 
stant magnification factors in the near-IR wavelengths 
(i.e., no differential magnification from 1.1/xm imaging 
data with HST/'WF CZ to S'pfteer/IRAC at 4.5/rm, see 
ICalanog et ahllMl^ . These factors, as well as the far-IR 
magnification factors, are shown in Table [5] The mea¬ 
sured flux densities and errors can be found with values 
from other wavelengths in Tabled 


4. MAGPHYS AND SPECTRAL ENERGY DISTRIBUTIONS 

The IRAC photometry measured in the 3.6 and 4.5/im 
bands, after carefully disentangling the foreground from 
background components, for the galaxies at z = 1 to 3, 
corresponds to rest-frame optical to near-infrared wave¬ 
lengths. These measurements redwards of the B aimer 
break allow us to measure the stellar mass and to break 
certain degeneracies in the SED modeling of DSFGs 
(though even with excellent photometry, the determina¬ 
tion of stellar masses is still degener ate with respect to 
assum ed star-formation histories, see iMichalowski et al.l 
I20I2I and IMichalowski et al.l 1^1411 . To model the spec¬ 
tral energy distributions of these DSFGs, we used the 
Multi-wavelength Analysis of Galaxy Physical Proper¬ 
ties (M agphys) program developed by Ida Cunha et al.l 
(|2008[1 . Magphys is a model package that empirically 
derives galaxy properties based on observations at rest 
wavelengths in the range 912A < A < 1mm. 

The SED models assume a Chabrier (2003) initial 
mass function (IMF) that is cutoff below 0.1 and above 
1OOM0. If a Salpeter IMF is assumed instead, we find 
that the stellar masses would be a factor of at least 1.5 
larger. The SED models do not include any AGN com¬ 
ponents, which may contribute significantly to the con¬ 
tinuum mid-IR emission. By ignoring AGNs, we could 
be overestimate the stellar mass of th e galaxies by as 
much as 0.3 dex (jda Cunha et al.l[2015[l . The templates 
are based on the Bruzal & Chariot (2003) spectrum li¬ 
brary used as an input for the SED models. These spec¬ 
tra assume an underlying continuous star-formation rate 
history described by a formation age and a time-scale 
parameter with random bursts superimposed on the con¬ 
tinuous model. We used flux density values from l/tm to 
880/(m shown in Table |3] as data to be modeled by Mag¬ 
phys. The resulting models selected based on Bayesian 
approach are able to determine the likelihood functions 
for optical depth of the dust, temperatures of the cold 
and warm dust components, stellar mass, star-formation 
rate averaged over the last 100 Myr, and dust mass and 
luminosity. 

5. RESULTS AND DISCUSSION 

A full summary of the fitted parameters for the six 
lensed DSFGs is found in TablelH We also show the best- 
fit total and unattenuated spectral energy distributions 
for each of the DSFGs in Figured 

Since rest-frame absolute H-band magnitudes (Mh) 
provide a guide to galaxy stellar masses that is not 
complicated by details of the assumed star-formation 
history, we compare the predicted values based on 
the best-fit SED models to other samples in Figure 
2. T he data comes from 850 and 870um-s elected 
SMGs (jHainline et al.ll20lil: iWardlow et al.ll2nilh . The 
Herschel-selecieA lensed galaxies are on average consis¬ 
tent with the rest-frame H-band magnitude distribu¬ 
tion of 850 and 870/rm-selected SMGs. Using a stan¬ 
dard Kolmogorov-Smirnov test on the two populations 
of SMGs identified from our Herschel observations and 
that of 850 and 870/tm observations, we find a K-S statis¬ 
tics of 0.253 and a probability value of 0.648, suggesting 
that the two data sets are from the same parent popula¬ 
tion (null hypothesis). Thus if there is no active galactic 
nucleus contribution to their H-band luminosities, and 
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Table 3 

Intrinsic Flux Densities of SMGs 



HBootes02 

HFLS08 

G09vl.40 

G12v2.30 

G15v2.19 

HLockl2 

i^Fiosw (z^Jy) 

— 

— 

— 

— 

1.63 ±0.9 

— 

^Fiiow (/^Jy) 

— 

0.7 ± 0.2 

— 

— 

2.6 ± 0 . 3 '' 

3.5 ± 0.7 

(MJy) 

— 

— 

— 

0.10 ±0.02^ 

3.4 ±0.5^ 

— 

-^Fieow (z^Jy) 

— 

— 

— 

— 

11.0 ± 1.1 

— 

(z^Jy) 

— 

— 

— 

— 

14.2 ±2.5 

— 

-Fks (MJy) 

2.5 ± 1.4 

— 

1.31 ±0.22 

0.74 ±0.05^ 

12.4 ± 1.7 

— 

^3.6Mm (ZiJy) 

9.0 ± 2.4 

12.2 ± 2.5 

3.7 ±0.3 

5.0 ±0.4 

77.9 ±6.2 

20.9 ± 2.1 

F 4 .5^111 (/^Jy) 

18.1 ± 4.8 

12.6 ± 2.6 

4.2 ±0.4 

9.0 ±0.8 

63.9 ±5.0 

37.7 ± 3.8 

^S250 (mJy)^ 

15.4 ± 1.0 

57.3 ± 6.67 

25.4 ±0.7 

38.0 ±1.3 

64.4 ±4.6^ 

149 ± 6.7 

Fs350 (mJy)^ 

18.9 ± 1.0 

62.0 ±6.67 

24.9 ±0.7 

46.8 ±1.3 

34.0 ± 2.7^ 

106 ±6.7 

^ssoo (mJy)^ 

15.1 ± 1.0 

44.7 ±6.67 

15.8 ±0.7 

38.8 ±1.3 

14.7 ± 1.o'’ 

52.7 ±6.7 

Fssso (mJy) 

4.11 ± 0.82^ 

— 

4.01 ± 0.94^ 

11.3 ± 2 . 3 '* 

— 

— 


Note. — De-magnified fluxes and photometric errors, which includes magnificat ion errors. The s e valu es were used as inputs for Magphys 
model fitting. Flux densities and errors in FllOW, H, and Ks bands are from fCalanog et al.l II2014I) unless otherwise noted. Redshift 
measurements were obtained using CO observations. 

Apparent flux densities and errors from iFu et al.l l|2012l') . 

Flux densities and errors from lMessias et al.l II2014I) . 

Ty pical errors in SPIRE photometry are about 10 mjy (magnified), which includes both statistical and confusion noise IlSmith et al.l 
12012). At far-IR waveleng ths we assume the foreground lenses are non-dusty galaxies, so the fluxes are due to the background sources only 
fsee IWardlow et al.ll201^ . The exception is G15v2.19, which is lensed by an edge-on, star-forming spiral galaxy. 

Flux densities and errors from IBussmann et al.l 1120131) . Values presented do not include absolute flux density calibration uncertainties of 
7%. 


Table 4 

Magphys Results of SMGs 


Name 

/^(SFH/IR)' 

TV 

2 ^ism 

pK] 

/7-1BC 

-^w 

[K] 

0 

0 

SFR 

[M0yr-l] 

sSFR 

[10-l°yr-l] 

Adust 

[10“Lo] 

A^dust 

[10«Mq] 

HBootes02 

Q 1-1-1-0.09 

Q o«.-|-0.12 

^•oJ-_0.69 

24.6ti;« 

41.7l?;° 

4± 1 

383li3 

9.3 ± 1.0 

'^•^-0.9 

15.2111 

HFLS 08 

n Q4+0 07 
G.oA _0 18 

^•'^^-0.13 

^•-L^_0.40 

24.9t3 4 

45.111.6 

0.9 ±0.1 

412+^*^ 

^l^_40 

49 ± 1.5 

13711° 

34.811" 

G09vl.40 

n qq-|-0.08 
U.OO_o.i6 

Q oy-l-0.16 
-0.16 

^•^^-0.40 


41.7111 

0.8 ±0.1 

1291°^ 

16.5 ± 1.0 

51.711,4 

8.99ll^ 

G12v2.30 

G.DD_o.o5 

q 7Q-I-0.26 

'^-0.03 

25 . 0 ^ 1:1 

53.111°, 

2.9 ±0.6 

lOllf 

50.7 ±0.6 

19011J 


G15v2.19 

0.68l°;?3 

0-53l°:“ 

^•'^^-0.19 

24.91°;" 

37.611“ 

1.8 ±0.2 

98.8lf 

5.3 ±0.5 

r)o 7-I-O.5 
' -0.4 

7 7cr-|-0.25 
'•'^-0.25 

HLockl2 

U.OO-O.OO 

^•^^-0.05 

5 10+0.84 

o-j-y_o.oo 

25.011° 

44.8111 

3.0 ±0.4 

495lf 

16.5 ±0.5 

202 IIJ 

29.711? 


Note. — is the fractional energy absorbed by the ISM, calculated from stellar-dominated (SFH) and dust-dominated (IR) photometry. 
TV is the total V-band optical depth of the dust seen by young stars in their birth clouds. is the equilibrium temperature of cold 

dust in the ambient ISM. is the equilibrium temperature of warm dust in stellar birth clouds. M* is the stellar mass. SFR is the 

star-formation rate averaged over 10® years, while ips = SFR/M* is the specific star-formation rate normalized to stellar mass. Adust 
Mdust fbe total dust luminosity and mass, respectively. Reported values are the best-fit values with their 16th and 84th percentiles. 

Values in the first and second rows give and respectively. 
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the star-formation histories are similar, they are likely 
to have stellar masses consistent with other starbursting 
galaxy samples. 

In Figure 3 we compare stellar mass, IR luminosity, 
and gas mass fraction of the sample to other Herschel- 
selected lensed gala xies from the literature, following 
iNegrello et al.l (1201411 . Among the six DSFGs presented 
here, we find gas mass values from CO molecular line 
observations f or two of the gala xies fG12v2.30 from 
iFii et al.1 120121 and G15v2.19 from iMessias et HI 120141 ). 
We make use of the stellar mass estimates from Mag- 
PHYS modeling to calculate /gas = Algas/(Mgas + M*) 
and find values of 0.20±0.06 and 0.20±0.07 for G12v2.30 
and G15v2.19, respectively. These are comparable to the 
gas mass fraction values quoted in the literature. 

We find that a significant attenuation by dust (ry Ri 4- 
5) is required to be consistent with the HST/WFC3, 
Keck/NIRC2, and IRAC photometry. Such high 
extinction values are consistent with other ULIRGs 


and SMGs (e.g. IGeach et al 

120071: Ida Gunha et alJ 

20101: IMichalowski et al.l 201C 

Hainline et al.l 120111: 

Wardlow et al.l 120111: IRowlands 

et al.l I2014I1. The stel- 


lar masses of these galaxies span the range of 8 x 10^° - 
4 X 10^^ Mq while their 100 Myr-averaged star-formation 
rates range from 100 to 500 Mq yr“^. These star- 
formation rates are consistent with the starburst nature 
of these DSFGs, some of which can also be classified 
as SMGs based on their lensing magnification-corrected 
sub-millimeter flux densities. 

The physical properties that we derive for our sample 
of lensed DSFGs clearly indicate fundamental differences 
between this population and other star-forming galaxies 
identified from optical surveys at similar or lower red- 
shifts. In Figure 0] we show the SFR (inferred from the 
SED) vs stellar masses for these galaxies and compare 
our sample of DSFGs to star-forming galaxies and SMGs 
from the literature. For reference we also plot the aver¬ 
age “main sequence” relations for galaxies at z = 1 and 
2. The horizontal lines show the selection, given that the 
magnification-corrected 500/rm flux density values have 
a minimum at the level of lOmJy for our sample. The 
selection suggests the redshift dependence of the star- 
formation rate vs stellar mass for our sample and other 
SMGs, since the selection based on a flux density results 
in selecting starbursts with higher star-formation rates 
with increasing redshift. 

Systems identified from the sub-millimeter surveys as 
SMGs have relatively large estimated stellar masses com¬ 
pared to the norma l star-forming galaxie s identified from 
standard surveys (|Magnelli et al.l [201 2tl . In fact, our 
lensed Herschel-selected DSFGs have SED-inferred stel¬ 
lar mass estimates in the range of ~ 8 x 10^° — 4 x 10^^ 
Mq with a mean value of ^ 2.3 x 10^^ Mq. This 
tight range of stellar masses is believed to not be re¬ 
lated to the optical and near-IR sensitivities for these 
systems, as lower mass galaxies are expect ed to be de¬ 
tected in the deep Spitzer observations (jReddv et al.l 
120061 : iMichalowski et al.l[2010tl . As mentioned earlier, we 
do not account for an AGN contribution to the SED 
of these galaxies. The current data do not allow us 
to separate AGN templates from SF in galaxies, while 
existing studies that separated AGN activity relied on 
significant mid-IR coverage, including Spitzer/IRS spec¬ 


tra (iBerta et al.l[2013[l . As discussed in Ida Cunha et all 
(l20m stellar mass is the single parameter that is most 
affected by an AGN contribution to the SED that is ig¬ 
nored in model fits, and we may have over-estimated 
the stellar mass by 0.3 dex. The star-formation rates 
measured for SMGs contribute very little to the stellar 
mass, with most of t he mass being built up before the 
sub-millimeter phase (IMichalowski et ^1201011 . This ex¬ 
plains the scatter that we see in the M*-SFR plane for 
the SMG population, with rates ranging from ^ 100 - 
400 Mq yr“^. 

While our lensed DSFGs have high stellar masses, com¬ 
parable to SMGs, their sSFRs may not amount to the 
high level seen in SMGs with values above 5 Gyr“^. 
Star-forming galaxies at z ^ 2 on the main sequence of 
that redshift have s SFRs at the level of 2 Gyr“^ or be¬ 
low ()Fu et al.l 1201311 . Based on SED models, only two 
of our targets have sSFRs that are at the level of 5 
Gyr“^ (G12v2.30 and HFLS08). It is yet unclear if all 
SMGs are a result of galaxy mergers, ho wever, a previ- 
ous multi-wavelength study of G12v2.30 (iFu et al.ll2012ll 
and the resulting lensing model did find clear evidence 
for a merger in the source plane of that system. Unfor¬ 
tunately, HFLS08 does not have the same level of multi- 
wavelength data as G 12v2.30 and the lens reconstruction 
(|Calanog et al.l2014fl for near-IR observations only shows 
a single galaxy. This is somewhat misleading since other 
lensed SMGs/DSFGs, including G12v2.30, show multi- 
component structure in the source plane with rest-frame 
optical emission spatially disjoint from the peak mm- 
wave dust em i ssion where the starburst is active (e.g., 
iFu et al.ll2012l : iDve et ani2015l : iHodee et al.ll201^ . 

The fact that only galaxies with sSFRs greater than 5 
Gyr“^ are SMGs does not imply that DSFGs with sS¬ 
FRs greater than 5 Gyr“^ are single galaxies. G15v2.19 
with a sSFR of roughly 1 Gyr“^ shows a clear merger 
with again distinc tively separate peak r e st-frame optical 
and d ust emission IMessias et al.l (1201411 ; iTimmons et al.l 
(|2015ll . It may be that a majority of the lensed DS¬ 
FGs are the result of merger-driven star formation activ¬ 
ities. Higher resolution observations are clearly desirable 
to address this, not just in the millimeter wavelengths 
with interferometers, but also in the optical and infrared. 
While Keck in the band and HST/WFC3 at 1.6^m 
and below can provide the necessary data to study DS¬ 
FGs in high resolution, our study finds that observations 
with Spitzer/IRAC are also desirable to determine key 
parameters from the SED, especially stellar mass and 
galaxy extinction. Such information can be extracted 
from the data despite issues related to the low spatial 
resolution in the data involving lensed systems. In the 
era of JWST we may be able to do more since imaging 
then at wavelengths around 3-4^m will provide necessary 
spatial resolution at the same level as interferometers. In 
that case we may not be able to just obtain globally av¬ 
eraged values, but also perform SED modeling of the 
clumpy interstellar medium of lensed DSFGs at spatial 
resolution of 300pc, as ac hieved by mm-wa ve interferom¬ 
eter data such as ALMA (jPve et al.ll2015f l and establish 
the stellar masses of individual clumps. 

6 . SUMMARY 

We have presented the rest-frame optical spectral en¬ 
ergy distribution and physical properties of six Herschel- 
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Figure 2. Magnification-corrected, rest-frame absolute H-h&nd 
magnitude Mh versus redshi ft for DSFGs. We also include 850 
and 8 70/im-selected galaxies JHainline et al.ll2QlU: IWardlow et al.l 
l20TiD. as well as Herschel-selected targets SDP.81 and SDP.130 
JHopwood et al1l201Ul . A Kolmogorov-Smirnov statistic of 0.253 
with a significance level of 0.648 suggests that Herschel-selected 
lensed galaxies are consistent with the rest-frame H-band magni¬ 
tude distribution of 850 and 870/im-selected SMGs. 


selected gravitationally lensed dusty, star-forming galax¬ 
ies (DSFGs) at redshifts of 1 to 3. The lensed DSFGs 
were first identified with Herschel/SPIRE imaging data 
from Herschel Astrophysical Terahertz Large Area Sur¬ 
vey (H-ATLAS) and Herschel Multi-tiered Extragalactic 
Survey (HerMES). 

1. The targets were observed with Spitzer/IHAC at 
3.6 and 4.5^m. Due to the spatial resolution of the 
IRAG observations at the level of 2", the leasing 
features of the background DSFG is blended with 
the flux from the foreground leasing galaxy in the 
IRAG imaging data. We make use of higher resolu¬ 
tion Hubble/WFCS or Keck/NIRG2 Adaptive Op¬ 
tics imaging data to fit light profiles of foreground 
leasing galaxies as a way to model the foreground 
components in order to disentangle the foreground 
lens and background source flux densities. The flux 
density measurements at 3.6 and 4.5/rm, once com¬ 
bined with Hubble/WFG3 and Keck 2.2/rm data, 
provide important constraints on the rest-frame op¬ 
tical spectral energy distribution of the Herschel- 
selected lensed DSFGs. We model the combined 
UV to mm-wavelength SEDs to establish the stel¬ 
lar mass, dust mass, star-formation rate, and vi¬ 
sual extinction, among other parameters for each 
of these Herschel-selected DSFGs. 
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Figures. Top: gas fraction /gas = Mgas/(Mgas -f Matellar) 
vs. redshift. We use gas mass values presented by |Fu et al.| 
1 I 20 TII and IMessias et al.( <120141) to determine gas fractions for 
G12v2.30 and G15v2.19, respecti vely. Gas fract i on val ue for the 
purple data point is provided bv IRiechers et ^ I|2Q1U) . Middle: 
magnification-corrected infrared luminosity Lir vs. redshift. Here 
we assume dust luminosities dominate the contribution to IR lu¬ 
minosities for our sources, as indicated by the blue points on the 
plot. Inf rared luminos i ty val ue for the purple data point is pro¬ 
vided bv IConlev et all 1120111) . Bottom: magnification-corrected 
stellar mass M* (bottom) vs. redshift. Inclu ded are SMGs—lense d 
and unlensed—exami ned in other studies: CoMev et al.l l|2011l) : 
IFu et ah 1 HOTD; and INegrello etHI ||2014I). Stella r mass 

value for the purple data point is provided bv Isolde v et al.l II2011I ') . 


2. Herschel-selected lensed galaxies are consistent 
with the rest-frame H-band magnitude distribution 
of 850 and 870/tm-selected SMGs. Assuming no 
contribution from AGN to their H-band luminosi¬ 
ties, these galaxies are likely to have stellar masses 
consistent with other starbursting galaxy samples. 
Herschel selections may also depend on redshift, 
since the selection based on flux densities results in 
selecting starbursts with higher SFRs with increas¬ 
ing redshift. 

3. The high extinction values for these DSFGs (ry ~ 
4-5) are consistent with other ULIRGs and SMGs. 
The stellar masses of these galaxies span the range 
of 8 X 10^° - 4 X 10^^ Mq, while their star-formation 
rates range from 100 to 500 Mq yr“^. The in¬ 
ferred SFRs are consistent with the starburst na¬ 
ture of DSFGs, some of which can also be classified 
as SMGs based on their flux densities. 

4. We find a large scatter between stellar mass and 
SFR for the SMG population. However, we also 
observe a correlation between the specific star- 
formation rate and dust temperature. This sug¬ 
gests that these galaxies were formed by merging 
systems. 

5. Lensed systems are intrinsically similar to far- 
infrared selected samples, which allows analyses of 
lensed galaxies to be generalized to the larger pop¬ 
ulation of SMGs and DSFGs, and vice-versa. 

6. We conclude that, despite typically low resolutions, 
Spitzer/IHAC data provide vital constraints on key 
parameters such as stellar mass. Followup observa¬ 
tions with higher resolutions by instruments such 
as the future JWST can not only introduce fur¬ 
ther constraints, but may also allow for detailed 
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Fig ure 4. Star-format i on rat e versus stellar mass o f SMGs and DSFGs. C ircles represent source s in this study, and we also include data 
from [Michalowski et al.l II2010I ]. IBanerii et al.l 1120111 ). IConlev et al.l (120111 ). IFu et all (120121 ). and IFu et aLl (120131 ). 850/im-selected SMGs 
are shown as diamonds and Herschel-selected targets are shown as stars. Red poi nts correspond to 1, blue points 2 : ~ 2, and green 
points z 3. The dashed lines show the main sequence correlations at each redshift dSueagle et al.|[2014l ). The horizontal dash-dotted lines 
represent the S 500 = 10 mjy selection cutoff for SMGs, drawn from mock SEDs at each redshift. 


analyses of the ISM of DSFGs with resolutions of 
^ 300pc, as achieved by previous interferometers 
such as ALMA. 
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7. APPENDIX 








Figure 5. Spectral energy distributions of each DSFG. Reported redshifts correspond to their respective sources. The red data points 
are the input fluxes (see Table [Sj. The blue and black curves are the best-fit attenuated and unattenuated SED’s, repectively. Redshifts 
correspond to the DSFGs and x? statistics indicate the goodness of fits for each SED. The lower panels display the residuals for each fit. 
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Figure 6. Spitzer IRAC de-blending results. Upper rows correspond to IRAC Channel 1 (3.6^m) and lower rows correspond to IRAC 
Channel 2 (4.5^m). From left to right we show: (a) high-resolution data, (b) the original Spitzer IRAC data, (c) the best-fit foreground 
model, and (d) the residual image with the foreground component(s) removed. We also show the xi statistic for each fit. Images are 
oriented such that North is up and East is left. Each cutout has size 20^' x 20'^ with each tickmark representing 1". 
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Figure [6] — continued. 






















